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Abstract 

The low-temperature specific heat of a superconductor M03SD7 with T c = 2.25 ±0.05 K has been 
measured in magnetic fields up to 5 T. In the normal state, the electronic specific heat coefficient 
7„,, and the Debye temperature Qd are found to be 34.5(2) mJ/molK 2 and 283(5) K, respectively. 
The enhanced 7 n value is interpreted due to a narrow Mo-4d band pinned at the Fermi level. The 
electronic specific heat in the superconducting state can be analyzed in terms a phenomenological 
two BCS-like gap model with the gap widths 2Ai/ksT c = 4.0 and 2A2/fcsT c = 2.5, and relative 
weights of the mole electronic heat coefficients 71/7,1 = 0.7 and 72/7,1 = 0.3. Some characteristic 
thermodynamic parameters for the studied superconductor, like the specific heat jump at T c , 
AC(T c )/7 n T c , the electron-phonon coupling constant, X e - P h, the upper H C 2 and thermodynamic 
critical H c q fields, the penetration depth A, coherence length £ and the Ginzburg-Landau parameter 
k are evaluated. The estimated values of parameters like 2Ao/£;_bT c , AC(T c )/7 n T c , N(Ep), and 
\ e -ph suggest that Mo3Sb7 belongs to intermediate-coupling regime. The electronic band structure 
calculations indicate that the density of states near the Fermi level is formed mainly by the Mo-4d 
orbitals and there is no overlapping between the Mo-4d and Sb-sp orbitals. 
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I. INTRODUCTION 

Several years ago, Bukowski and coworkers reported the superconductivity in an inter- 
metallic M03SD7I1I]. The superconducting phase transition of this material has been deter- 
mined from the magnetization and electrical resistivity on a single crystal to be 2.08 K 111]. 

n □ 

Based on the magnetization [1], an d point-contact (PC) Andreev-reflection data |2j, the au- 
thors of Ref. Q, [2I concluded the BCS-type behaviour of the superconductor. However, in 
the recent work, Dmitriev et al.^ had suggested that MosSby is not a trivial BCS (s-wave) 
superconductor but rather it has a (s+g)-wave or another unconventional pairing symmetry. 
More recently, Candolfi et al. [4J, have reported magnetic susceptibility, specific heat, and 
electrical resistivity down to 0.6 K and have interpreted their data within the frameworks of 
the spin fluctuation theory. Thus, the issue of pairing mechanism becomes debated and the 
proper superconductivity nature of Mo3Sb7 deserves studies in more details. Nonetheless, 
information about the superconducting properties of MoaSby is still scarce. Except for the 
superconducting transition, other parameters like upper H C 2 and thermodynamic critical H c q 
fields, the penetration depth A^, coherence length £ and the Ginzburg-Landau parameter k 
are not known yet. In order to determine these fundamental thermodynamic parameters of 
the MosSby superconductor, we have performed specific heat C v at temperatures 0.4 -10 K 
and in magnetic fields up to 5 T. In this work, we report the results of the C v measurements 
and we try also to address the question: which electron pairing model is more suitable for 
description of the superconductivity in MosSby? The specific heat data are supplemented 
by electronic band structure calculations performed using the the full-potential local-orbit 
method. 



II. EXPERIMENTS 

A powdered sample of MosSby was prepared from mixture of powder metals Mo and Sb 
with purities of at least 99.95% from Alfa Aesar by annealing in an evacuated silica tube at 
700° for three days. The mixture was then grounded and pressed into the form of cylinder, 
and annealed at 600° for two days. The purity, homogeneity and composition of the sample 
were checked by X-ray diffraction and dispersive spectroscopy. 

Fig. CTJ shows the experimental powder X-ray diffraction patterns of MosSby and as well 
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as the profile fit and difference patterns in the Rietveld analysis. It can be seen, the observed 
Bragg reflections were indexed on the base of the cubic h^Gey-type structure (space group 
Im-3m) 5j. In the refinements, we assumed that Mo atoms occupy the 12e (x = 0.342) and 
the antimony distribute on the 12d and 16f (x = 0.162) sites, respectively. The refinements 
of the X-ray intensities yield lattice parameter a = 0.9551 (1) nm, in well agreement with 
that previously reported by Dashjav et al. for a single crystal (0.9559 nm) [6j. 

The microanalysis was performed by collecting several EDX spectra at various locations 
across the surface of a sintered pellet. The compositions determined by EDX Mo: 29.2 
± 3 at.% and Sb: 70.8 ± 3 at.% are close to the ideal composition Mo:Sb= 3:7. No impurity 
phase within the experimental error is observed in the EDX spectra. The SEM image of the 
typical samples is shown in the inset of Fig. [TJ 

The specific heat was measured using a thermal relaxation method in the temperature 
range 0.4 - 10 K and in magnetic fields up to 5 T, utilizing PPMS system (Quantum Design). 
The experimental errors were taken as found by Marriott et al. 7J], i.e., to be less than ± 
2 % for T> 1.5 K and about ± 20% at 0.4 K. Before the specific heat study, the bulk 
superconductivity in the investigated sample was examined by means of magnetization and 
electrical resistivity measurements. From the latter measurement we obtained the residual 
resistivity in the normal state p n = 95x1 -8 fhn. 

The electronic band structure calculations were carried out using the full-potential local- 
orbit method (FPLO-5) 0,3,0. The calculations were performed for the reciprocal space 
mesh containing 816 points within the irreducible wedge of the Brillouin zone. 

III. RESULTS 

The calculated densities of states (DOS) of MosSby near the Fermi level are shown in 
Fig. El together with the partial contributions from the constituent Mo and Sb atoms. 



' Self-consistent 
6]. Both calcu- 



The global feature of the DOS is in reasonable agreement with the results o 
LMTO (linear muffin tin orbitals) calculations performed by Dashjav et al 
lations have indicated the dominating contribution of the Mo-4d orbitals to the total DOS 
The calculated DOS at the Fermi level N(Ep) is estimated to be 42 states/eV unit cell, cor- 



responding to about 1.04 states/eV atom. Comparison of the o 
those of YBa2Cus07 and MgB2 (~ 0.12 states/eVatom.spin 11]- 



)tained N(Ep) value with 



12| ) points to an enhance- 
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ment of the N(Ep) of M03SD7. This observation implies that the 4d-electrons of the Mo 
atoms forming a narrow band at the Fermi level give rise to the density of states, in similar 
manner as in the case of the A15 snpercondnctors Q. A distinct difference between onr 
and the previous published data lies in the fact that the position of the Fermi level in our 
calculations is shifted toward higher energy by 0.02 eV. Thus, Ep situates in the opposite 
site of a local valley. The important result emerging from our calculations is not only the 
fact that the largest peak of DOS located at -0.035 eV below Ep (band width Eb of ~ 0.05 
eV) forms mainly from the Mo-4d orbitals but also no overlapping between Mo-4d and Sb-sp 
orbitals is observed. Thus, the hybridization of the Mo-4d with the valence sp-band of the 
Sb atoms seems to be negligible and the Mo-4d and Sb-sp orbitals contribute to the total 
DOS probably separately. We recognize also that the contribution of the sp-band of the Sb2 
atoms to the total DOS at Ep is much lower than that of Sbl. This finding may support 
a very weak hybridization strength between Mo-4d and Sb-sp orbitals, since the distance 
Mo-Sb2 (~ 2.78 A) is shorter than the Mo-Sbl (~ 2.82 A) one. 

The specific heat C p data of MosSby at zero field is shown in Fig. [3] in the form C p /T 
vs. T 2 . The occurrence of superconductivity is evidenced by a specific heat jump at 2.2 
±0.05 K. The superconducting phase transition temperature determined from the specific 



heat measurement is consistent with those obtained from the magnetization and e 



resistivity (not shown here) and with the specific data reported by Candolfi et al. {4]. Note, 
the value of T c observed in the present sample is of about 0.1 K higher than that found in a 
single crystal of 2.08 K, jlj, though specific heat data for single crystal are not yet available. 
The difference in the superconducting transition values between the single-crystalline and 
polycrystalline samples seems to be associated with the change in the densities of states at 
the Fermi level. Since the unit cell volume of the studied sample with 871.34 A 3 is smaller 
than that of the single-crystalline sample (879.769 A 3 ) [l|, one supposes that the decrease 
in the volume results in displacement of the 4d- level to a higher energy respect to Ep, and 
in consequence increases DOS. 

The normal-state specific heat data of Mo3Sb 7 can be analyzed using the standard for- 
mula: 

C p /T = ln + (3T 2 , (1) 

where the first term gives the electronic contribution and the second represents the lattice 
contribution to the total specific heat. For the temperature range 2.5 - 10 K, it is found 



ectrical 



that 7„ = 34.5(2) mJ/molK 2 and (3 = 0.85(5) mJ/molK 4 . The latter value corresponds to 
the Debye temperature of 283(5) K. 

From the specific heat data (Fig. [3]), the dimensionless specific heat jump AC(T)/7 n T c 
was evaluated to be 1.56. Thus, there is a substantial difference in the AC(T)/^ n T c values 
obtained by us and by Candolfi et al. 4] i.e., 1.04. Since, the value of 1.43 is the weak- 
coupling limit [14J], the superconductivity in the sample studied here could belong to the 
strong-coupling regime. An information about the strength of the electron-phonon coupling 
can be gained by evaluating the average electron-phonon coupling constant X e - P h from the 
McMillan equation 151 ]: 

= 1.04 + /x*Zn(8 D /1.45T c ) 

(1 - 0.62/z*)/n(e D /1.45T c ) - 1.04 [ ' 

where fi* is the Coulomb pseudopotential and it is usually taken between 0.1 and 0.15 
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16|. Putting n* = 0.15, Q D = 283 K and T c = 2.2 K into Eq. [2] we obtained A 



e—ph 



0.57. For MosSb7 the value of fi* ~ 0.17 can be roughly estimated from the electronic band 
gap E B = 0.05 eV, phonon frequenc y a;n ~ &d = 283 K and product N(E F )V ~ 0.2 using 



the equation: fx* = TT^rw^njE I — ) — 1- ^he value of the product N(E F )V was calculated 



l+N(E F )Vln(E B /uj ) 

with help of the BCS equation: T c = l.UQ D exp(-l/N(E F )V). 

The knowledge of 7„ allows us to estimate the electron-phonon mass enhancement A 7 
parameter, which should be similar in magnitude to \ e - P h- Employing the relation 7 = 
|7r 2 iV J 4/c^(l + \.y)N(Ep) and taking N(E F ) as the DOS obtained by the electronic band 
structure calculations, we obtained A 7 = 0.64. Compared to weak-coupling \ e - p h value of 
0.4, the relatively large and consistent values of X e -ph an d A 7 suggest an intermediate strong 
coupling strength. This finding is consistent with an enhanced value of the specific heat 
jump at T c . 

We derived the electronic specific heat C e \ = C p — /3T 3 and have plotted the data in Fig. 
|4] a in the form of semi-logarithmic scale ln(C e i) vs 1/T. Clearly, in the superconducting 
state, the C e i(T) curve reflects an exponential behaviour, indicating a gap character of the 
studied superconductor. We wish to add that within the experimental error we have not 
seen any power law T 2 or T 3 in the temperature dependence of the electronic specific heat. 
This observation rules out any anisotropic symmetry of superconductivity as suggested by 
Dmitriev et al. . However, an attempt to fit the experimental data with the phenomeno- 
logical relation for the conventional s-wave superconductors: C e \ = %.b^ n T c exp{— 1.44T C /T) 
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[17(|, for the whole temperature range studied (represented as dashed line) is failed. In fact, 
the single-gapped BCS-like C e i = A^ n T c exp{— A //cbT) function describes well the exper- 
imental data only in a limited temperature range, e.g, for 1 K < T < 2.2 K we obtained 
parameters A = 12.7 and A /k B = 4.0, corresponding to 2A //cbT c 3.63. 

In order to reproduce the data in the whole temperature range studied, one applies a 
phe_.ogica. two-gap model, which was proposed by Wang et al. Q, and by Bouquet 
et al. [19( some years ago for MgB2. Following the latter authors, the electronic specific 
heat below T c can be considered as the sum of the contributions of two bands characterized 
by respective gap widths Ai and A2, and by partial Sommerfeld ratio 71 and 72, where 
7n, = 71 + 72- We have fitted the experimental data to theoretical curves within the a-model 



20j . The theoretical data were calculated based on the following equations: 

C el {t)hT c = tid/dt)^)/^) (3) 

roc 

S el {t)hT c = -(3a/7r 2 ) / ]/(x)ln(/(x)) + (1 - /(x))ln(l - f(x))]dx (4) 







where a is the gap ratio a = A(T)/ Abcs(T) with 2A j bc*s , (0) = 3.5ksT c and f(x) 



[exp( a ^ +s ® 1 ) + 1] x . 5(t) is the nornm 



t = T/T c and tabulated by Muhlschlegel 2l|. The best fitting result was observed for 



ized BCS gap at the reduced temperature 



the following parameters 2Ai = 4.0A;bT c and 2A2 =2.5A;bT c and relative mole weight 71/72 
= 70/30. The good fitting of the specific heat data to the exponential law implies that 
Mo3Sb7 may be an s-wave-like superconductor with two energy gaps dependent on temper- 
ature. The existence of two bands having different electronic densities at the Fermi level can 
be understood if we assume that a narrow Mo-4d band and a broader Sb-sp band contribute 
differently to N(Ep). Based on the partial gap widths, one suspects that the larger gap 
determined by the stronger electron-phonon coupling would correspond to the Mo-4d band, 
while the smaller gap associates with the Sb-sp band, and then with 7 n = 34.5 mJ/molK 2 
and the relative 71/72 ratio, one estimates the partial electronic coefficient 72 =1.5 mJ/(g 
at.Sb K 2 ) and 71 = 8.0 mJ/(g at. Mo K 2 ). This viewpoint is quite accordable with the results 
of the electronic band structure calculations shown above, where the dominating input of 
the Mo atoms to the DOS is observed. We may add that two-band gap feature has been 



found in some superconductors, except for MgB2 mentioned above, also in YbNi2B2C 22] 



and Nb^Sn 



23|]. 



The electronic entropy in the superconducting state was calculated through Ssc 



Jq C C e i/TdT (Fig. H]b). At T c , S$c amounts to 80 mJ/molK and it is comparable with 
the entropy in the normal state S n = 7 n T c of 76 mJ/molK assuming a constant 7 n = 34.5 
mJ/molK 2 below T c = 2.2 K. 

The temperature dependence of the thermodynamic critical field HoH c0 is obtained by an 
integration of the data in the superconducting state according to the equation: 

l -^V M Hl(T) = AU(T) - TAS(T) (5) 

where Vm is the volume per mole, AU is the internal energy difference and AS" is the entropy 
difference. From the equations: 

AU(T) = £\c el (T') - ln T')dT' (6) 

AS{T) = f C C ei( T ')-lnT' dT , (?) 
Jt t 

we computed the temperature dependence of AU and TAS (see inset of Fig. [5]). The critical 
field HqH c o is calculated from Eq. [5] is plotted in Fig. [5] as a function of temperature. An 
extrapolation to K yields fi H c0 (0) = 0.029 T. A theoretical (j, H c0 (0)- value can be approx- 



imated from the relation: ^ H c0 (0) = 7.65xl0" 4 7 1 / 2 T c [24(, to be w 0.027 T. This value is 
in good agreement with the experimental one. Using the relationship between the /io-ff c o(0) 
and the superconducting energy gap A(0): HqVmH^q (0) = (■^^ T )A(0) 2 we deduce A(0) to 
be 0.39 meV (~ 4.18 K). From this value, the gap 2A(0)//cbT c for MosSb7 is calculated to 
be 3.8, slightly larger than BCS value of 3.52 for weak-coupling superconductors. 

In order to analyze the upper critical field H c2 , we measured specific heat at several 
selected magnetic field strengths (Fig. E]). With increasing magnetic fields up to 2 T, the 
superconducting transition shifts down to lower temperatures. The resulting upper critical 
field against temperature is shown in Fig. [7] a. Apparently, the values of fioH c2 (T) can 
be described by a relation fi H c2 (T) = fi Q H c2 (0)[l - (T/T C ) L29 ] with fi H c2 (0) = 2.24 T 
and T c = 2.2 K. Near T c the initial gradient —fiodH c2 (T)/dT amounts to 1.25 T/K. An 
independent determination of —fi dH c2 (T)/dT may be made using the equation given by 



Orlando et al. 25j: — fiodH c2 (T) / dT = 4.48xl0 3 7po i n the dirty limit. For the experimental 
values of 7 and po we estimated —fiodH c2 (T)/dT w 1.12 T. Thus, there is satisfactory 
agreement with the experimental value. 
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Using the relationship poH C 2(0) ~ 0.69(— podH c2 (T) / dT)T c for a type-II superconductor 



in a dirty limit [26j, we estimated poH C 2(0) to be 1.9 T which is, however, slightly smaller 
than the value PqH C 2(0) fitted above. A better agreement between the theoretical and exper- 
imental values might have been obtained if one takes an additional parameter characterizing 
spin-orbital coupling X so . 

The value poH C 2(0) can be estimated with the help of the Maki theory j^]. The upper 
critical field is given through the relation: p H c2 (0) = aMH po /\/2, where the BCS value 
of the Pauli limiting field is H^ 3 = 1.84T C and olm is the Maki parameter. Taking into 
account the electron-phonon corrections to_the Pauli limiting process the field H po can be 



written as H po = H^ 3 ^1 + X e - P h 25|, [28|. The parameter olm is obtained through the 
relation % = -0. 528 fiodH^/dT or alternatively olm = 3e 2 fry n p n /(2m7r 2 /c^), where m is the 
mass of the free electrons. With p dH c2 /dT '=-1.25, and 7„ = 263 J/m 3 K 2 and p n = 95xl0 -8 
ffrn, the value of the Maki parameter amounts to 0.66 and 0.59, respectively. Using the 
average value olm = 0.625 and \ e - p h =0.6 we obtain the upper critical field fioH c2 (0) = 2.26 
T, which is close to the experimental value of 2.24 T. 

Having upper and thermodynamic critical fields, we can estimate the type-II characteristic 
Ginzburg-Landau (GL) parameters at K. From the relation: fi H c2 (0) = 2?r ^° , where $o 
= 2.067xl0 -15 Wb is the flux quantum we obtained coherence length ^gl ~ 12 nm. Using 
the relation: fi H c2 (0) = ^/2k{i H c (Q) we get Ginzburg-Landau parameter k = 54.6. The 
relation: fi H c (0) = 27i y| A 2 results in the penetration depth A ~ 660 nm. These parameters 



can be evaluated by other different method with the help of relations given in Ref. [25]. 
From A = 6.42xl0~ 4 p/^c and k = 2.37xl0 6 7 1 / 2 p, the parameters n and A were calculated 
to be 420 nm and 36.5, respectively. We suppose that smaller values of these parameters 
deduced from the resistivity may be due to an underestimated value of p. 

The field dependence of the Sommerfeld coefficient at 0.4 K, ^qak{.H) = C p /0A(H), 
is displayed Fig. [7] b. We observe linear dependence of the 70.4^ with p H up to 1.5 T, 
shown by the dashed line. The linearity of 7 .4^(i?) curve of Mo3Sb 7 accords with the 



theory for gapped superconductors 



29] . Although, such a linear behaviour of ^q A k{H) is 



not observed in fields up to PqH c2 . Above 2 T the joAxiH) levels off to a constant. Thus, 
the spin-fluctuation scenario in Mo3Sb 7 is not confirmed by our specific heat data, where two 
characteristic features of spin fluctuations, namely C P (T) ~ T 3 ln(T) 30j, and ^{H) ~ H 2 



3l|, are not found. 
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TABLE I: Characteristic parameters of the superconductor M03SD7. 



T c (K) 


2.2(0.05) 


7n [mJ/(molK 2 )l 


34.5(2) 


8 [mJ/(molK 4 )l 

/ L / V /J 


0.85(5) 


®D (K) 


283(5) 




0.6(0.05) 


AT(£p) [states/ (eV at.)] 


1.0 


AC(T c )/( 7n T c ) 


1.56 


2Ai/(A; B T c ) 


4.0 


2A 2 /(fe jB r c ) 


2.5 


-Un(dH r2 /dT) (T/K) 


1.25 


LLnHrO (T) 


2.24 


/xoi^co (T) 


0.029 


^GL (nm) 


12 


At 


36 - 55 


Al (nm) 


420 - 660 



IV. CONCLUSIONS 

In the present work we have presented a detailed analysis of the low-temperature specific 
heat data of the superconductor M03SD7. The lattice specific heat is characterized by the 
Debye temperature of 283 K. The observed Sommerfeld coefficient of 7 n = 34.5 mJ/molK 2 
and a large electronic density of states, N(Ep) of 1.04 states/eV atom, are indicative of 
a narrow band of the Mo atoms at the Fermi level. We have observed that the electronic 
specific heat in the superconducting state follows a single gapped function in a narrow 
temperature range. The gap width of 2A //c_bT c 3.64 is inferred for the data above 1 K. 
Assuming the presence of two superconducting BCS-like temperature dependent gaps we 
are able to describe the experimental data in the whole temperature range studied. The 
gap widths 2Ai/fcfiT c = 4.0 and 2A 2 / ksT c = 2.5, and relative mole weights of the mole 
electronic heat coefficients ji/j n — 0.7 and 72/ /r y n = 0.3 were fitted. From the analysis of the 
specific heat data we obtained characteristic thermodynamic parameters, like the specific 
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heat jump at T c , AC(T c )/7 n T c , electron-phonon coupling constant, X e - P h, upper H c2 and 
thermodynamic critical H c q fields, Ginzburg-Landau penetration depth Al, coherence length 
£ and Ginzburg-Landau parameter k. The obtained value of X e ~ P h, N(Ep), AC(T c )/7 n T c 
and 2A //cbT c can be accounted for intermediate electron-phonon coupling. We found a 
linear field dependence of the Sommerfeld coefficient at 0.4 K for fields below ~ 2 T. The 
data of electronic band structure calculations show no overlapping between Mo-4d and Sb- 
sp orbitals and therefore, may support the proposed two-gap feature. In order to confirm 
the two-gap phenomenon in Mo 3 Sb 7 , further studies using sophisticated techniques e.g. 
by means of spectroscopy techniques such as tunneling-, point-contact- or photoemission- 
spectroscopy and muon spin relaxation are planned. 
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FIG. 1: Powder X-ray diffraction pattern of M03SD7. The observed (open circles) and calculated 
(solid line) profiles are shown on the top. The vertical marks in the middle are calculated positions 
of Bragg peaks. The line in bottom of the plot is the difference between calculated and observed 
intensities. The SEM image of the typical samples is shown in the inset. 

FIG. 2: (Color Online) Total (dotted line) and partial densities (Sbl: dash-dotted, Sb2: dashed 
and Mo: solid line) of states of MosSby. 

FIG. 4: (Color Online) a) the electronic specific heat in a logarithmic scale as a function of 1/T. 
The dotted and dashed lines are theoretical ones (see text) The error bars are given with values 
of ± 20% at 0.4 K and ± 2% at 1.7 K. b) Temperature dependence of the electronic entropy at 
and 2 T. 

FIG. 5: (Color Online) The thermodynamic critical field H c q(T) as a function of temperature. 
The dotted is a guide to the eye. The inset shows temperature dependence of the internal energy 
difference AU and entropy difference multiplied by the temperature TAS. The dotted line is a 
guide to the eye. 

FIG. 6: (Color Online) The temperature dependence of the specific heat divided by temperature 
Cp/T at several magnetic fields. The dotted line is a guide to the eye. 

FIG. 7: (Color Online) a) Upper critical field H C 2 vs. temperature. The dashed line represents 
the fit of the experimental data to noH&^T) = ^0-^2(0) [1 — (T/T c ) 1,29 ]. b) Field dependence of 
the Sommerfeld coefficient at 0.4 K. The dashed line illustrates a straight linearity of the 'Joak(H) 
dependence. 

FIG. 3: (Color Online) The specific heat of MosSby divided by temperature as a function of square 
temperature. The dotted line is a guide to the eye. The dashed line is a fit to Eq. [TJ 
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